arXiv:1507.04143vl [stat.AP] 15 Jul 2015 


A Shock Model Based Approach to Network Reliability 


S. Zarezadeh* S. Ashrafi^^ and M. Asadi^ 


Abstract 

We consider a network consisting of n components (links or nodes) and assume that the 
network has two states, up and down. We further suppose that the network is subject to shocks 
that appear according to a counting process and that each shock may lead to the component 
failures. Under some assumptions on the shock occurrences, we present a new variant of the 
notion of signature which we call it t-signature. Then t-signature based mixture representations 
for the reliability function of the network are obtained. Several stochastic properties of the net¬ 
work lifetime are investigated. In particular, under the assumption that the number of failures 
at each shock follows a binomial distribution and the process of shocks is non-homogeneous 
Poisson process, explicit form of the network reliability is derived and its aging properties are 
explored. Several examples are also provided. 

Keywords: signature, fatal shocks, counting process, nonhomogeneous Poisson process, two- 
state networks, stochastic ordering. 


1 Introduction 

Networks include a wide variety of real-life systems in communication, industry, software engineer¬ 
ing, etc. A network is defined to be a collection of nodes (vertices) and links (edges) in which some 
particular nodes are called terminals. For instance, nodes can be considered as road intersections, 
telecommunications switches, servers, and computers; and examples of links can be telecommuni¬ 
cation fiber, railways, copper cable, wireless channels, etc. 

According to the existing literature, a network can be modeled by the triplet N = (V,E,T), in 
which V shows the node set, where we assume \V\ = m, E stands for link set, with \E\ = n, and 
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T C y is a set of all terminals. When all terminals of the network are connected to each other, 
the network is called T—connected. We assume that the components (links or nodes) of a network 
are subject to failure, where the failure of the components may occur according to a stochastic 
mechanism. A link failure means that the link is obliterated and a node failure means that all links 
incident to that node are erased. Assuming that the network has two states up, and down, the 
failure of the components may result in the change of the state of the network. 

In reliability engineering literature, several approaches are proposed to assess the reliability of 
a network. An approach, to study the reliability of a network with n components, is based on 
the assumption that the components of the network have statistically independent and identically 
distributed (i.i.d.) lifetimes Ai, A 2 ,..., A„, and the network has a lifetime T which is a function 
of Xi,... ,Xn. An important concept in this approach is the notion of signature that is presented 
in the following definition; see [IT] and |2|. 

Definition 1 . Assume that tt = (ej^, 6 * 2 , • • •, is a permutation of the network components 
numbers. Suppose that all components in this permutation are up. We move along the permutation, 
from left to right, and turn the state of each component from up to down state. Under the 
assumption that all permutations are equally likely, the signature vector of the network is defined 
as s = (si,..., Sn) where 

Si = niln\, i = l,...,n 

where Ui is the number of permutations in which the failure of ith component cause the state of 
the network changes to a down state. In other words, s* is the probability that the lifetime of the 
network equals the ith ordered lifetimes among Aj’s, i.e., Si = P{T = Xi,n), where Xi,n is the ith 
order statistic among the random variables Xi,X 2 , ..., X^. 

The signature vector depends on both the structure of the network and how to define its states. 
However, it does not depend on the real random mechanism of the component failures. Under this 
setting, the reliability of the network lifetime T, at time t > 0, can be represented as 

n 

P{T>t) = Y,SiP{Xi,n>t), (1) 

i=l 

see [18]. In recent years, a large number of research works are reported in the literature investigating 
different properties of the reliability function ([T|). We refer, among others, to |18j-[24] and references 
therein. 

Another approach, in assessing the reliability of a network, is recently proposed by Gertsbakh 
and Shpungin [9]. These authors consider a network with n components, and assume that the 
component failures appear according to a renewal process {N[t),t > 0} defined as a sequence of 
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i.i.d. non-negative random variables (r.v.s) Yi,Y 2 , ■■■,Yk ,.... The random variable N{t) shows the 
number of components that fail in the network on interval [0,t], and the failures in {N(t),t > 0} 
appear at the instants Sk = ^ — 1,2,.... Under the assumption that all orders of 

component failures are equally likely, the reliability function of the network lifetime T can be 
represented as 

n 

P{T >t) = Y^ S,P{N{t) =i), t> 0, (2) 

i=l 

where Si = ^k-,- Motivated by this, under the assumption that the failure of the network 

components occur according to a counting process, Zarezadeh and Asadi [22] investigated various 
properties of the model in ([2|) based on different scenarios. Zarezadeh et al. [23| studied stochas¬ 
tic properties of dynamic reliability of networks under the assumption that the components fail 
according to a nonhomogeneous Poisson process (NHPP). 

The aim of the present study is to give new models for the reliability of the network under the 
assumption that the components of the network are subject to shocks. We consider a two-state 
network and assume that the network is subject to shocks that appear according to a counting 
process. We further assume that each shock may lead to component failure and consequently the 
network finally fails by one of the arriving shocks. The reset of the paper is organized as follows: 
In Section 2, we obtain the mixture representations for the reliability of the network lifetime. For 
this purpose, a new variant of the notion of signature, call it t-signature, is introduced which allows 
us to assume that at same time more that one component failure may occur. We then compare 
the t-signature based reliability of two different networks under various assumptions. In Section 
3, we assume that the number of failed components in each shock are conditionally distributed as 
binomial distribution. Under this condition, mixture representations for the reliability function of 
the network are obtained and stochastic and aging properties of the network lifetime are investi¬ 
gated. In particular, we show that when the shocks arrive according to a non-homogeneous Poisson 
process (NHPP) and the arrival time of the first shock has increasing hazard rate average (IHRA), 
then the distribution of the network lifetime is IHRA. Section 4 is devoted to the reliability of 
the network under fatal shocks. It is assumed that at time of occurrence of a shock at least one 
component of the network fails. Under this assumption a mixture representation for the network 
reliability is obtained based on a new variant of the notion of signature. 

2 Network reliability under shock models 

In this section, we assume that the network is subject to shocks that appear according to a counting 
process. In reality, this may happen as a result of a sequence of heavy road accidents, floods. 
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earthquakes, fires etc. We explore the reliability of the network where each shock may lead to 


the failure of the network components. Before doing so, we define a variant of the concept of 


signature which avoids the restriction of not allowing the ties. To be more precise, let Xi, ...,Xn be 
i.i.d random variables representing the component lifetimes of the network. One of the assumptions 
that is necessary to define the notion of signature is that there do not exist ties between Xi ,..., Xn, 
i.e. P{Xi = Xj) = 0 for every i ^ j (see, for example, |20]1. However, in real life situation, this is 
possible that more than one component may fail at each time instant, i.e. ties may exist between 
Xi,... ,Xn- For example when the network is under shock, each shock may results the failure 
of more that one component at the same time. Under this assumption, in the sequel, we define 
a variant of the notion of signature. First let us define the discrete random variable M as the 
minimum number of components that their failures cause the network failure. Obviously M takes 
values on {1, 2,... , n}. Suppose further that n* is the number of ways that the components fail 
in the network and n* is the number of ways of the order of component failures in which M = i. 
Assuming that all the number of ways of the order of component failures are equally likely, we 
define the ’’tie signature” {t-signature) vector associated to the network as ..., where 

r rii . 

Si = —:> * = 

n* 

It should be noted that t-signature, similar to the concept of signature, depends only on the 
structure of the network and does not depend on the random mechanism of the component failures. 

In the following example, we compute the t-signature vector for a simple network. 

Example 1. Consider a network with 3 links and 3 nodes depicted in Figure [TJ The links are 
subjected to failure and nodes a and c are considered as terminals. We assume that the network is 
functioning if and only if terminals are connected. 
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Figure 1: Network with 3 links and 3 nodes. 


Let TT denote the order of link failures in the network . All possible vr and the associated M are 
presented in Table [H where the numbers in the braces indicate that the corresponding links failed 
at the same time. Hence, n* = 13 and the elements of the t-signature are calculated as 
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Table 1: All ways of order of links failures 


TT 

M 

TT 

M 

TT 

M 

(1,2,3) 

1 

({1,3},2) 

1 

({1,2,3}) 

1 

(1,3,2) 

1 

({2,3},!) 

2 



(2,1,3) 

2 

({1,2},3) 

1 



(2,3,1) 

2 

(3,{1,2}) 

2 



(3,1,2) 

2 

(2,{1,3}) 

2 



(3,2,1) 

2 

(1,{2,3}) 

1 




It is interesting to note that the signature vector of this network equals s = (i, |,0). 

The following lemma gives a formula for computing n*. 

Lemma 1. Let a n-component network he under shocks. Let n* be the number of ways that the 
components the network fail under the assumption of ties. Then 

n j ^ .X 

"• = EEu 

j=l k=0 ^ ' 

Proof. We use the following combinatorial argument; The number of ways to put n distinct objects 
into m distinct boxes, n > m, such that every box contains at least one object is 

k=0 ^ ' 

Let J be the number of shocks such that in occurrence of each one at least one component fails. It 
is clear that J takes value on {1,n}. If J = j, is fixed, the number of ways that the components 
numbers {1,2, ...,n} can be under j shocks is the same as the number of ways to put n distinct 
objects into j distinct boxes such that every box contains at least one object. Thus, summing up 
over j, j = I,..., n, we get 

n j y 

"• = EEu 

j=l k=0 ^ ' 

□ 

Consider a two-state network with lifetime T which is subject to shocks, where shocks appear 
according to a counting process, denoted by {f{t),t > 0}, at random time instants -di,'d 2 , • • • • 
We assume that each shock may lead to component failures and further assume that the network 
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finally fails by one of these shocks. Let random variable Wi, i = 1,... ,n, denote the number of 
components that fail at the fth shock and Wq = 0. If (t) denotes the total number of components 
that fail up to time t, then N{t) takes values on {1,2,..., n} and 

m 

N{t) = Y,Wi. 

i=0 

Under the assumption that the process of occurrence of the shocks is independent of the number 
of failed components, using the law of total probability, the distribution function of N{t) can be 
written as 

OO 

P{N{t) < x) =^P{N{t) < xm = k)Pm = k) 

k=0 

OO C(t) 

=E^(E W, < x|e(t) = k)P(^(t) = k) 

k=0 ^=0 

OO 

= Y^Hk{x)P{m = k), (3) 

k=0 

where Hk{x) denotes the distribution function of r.v. these assumptions, the network 

fails if N{t) > M. Hence, the network lifetime can be defined as 

T = min|A^(t) > Mj 
t>o 

and thus, we have P{T > t) = P{N{t) < M). Therefore, using the law of total probability and the 
fact that the total number of components that fail up to time t is independent of the t-signature, 
we get 


P{T > t) =P{N{t) < M) 

n 

= P{M = i)P{N{t) <i-l) 

i=l 

n 

= Ys'lP{N{t)<i-l). (4) 

i=l 

Let SJ = Yl'l=j+i ^1) then using ([3|) and (jH, we have 

n OO 

p{T >t)=Y^JYl 

i=\ k=0 

OO 

k=0 
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where, for fc = 0,1,... 


(3k,n = '^slHk{i - 1) 

i=l 

n—1 k 

=E'5jnE^* = •?■)• (6) 

j=0 i=0 

In the following proposition some properties of 13k,n are investigated. 

Proposition 1. Let be the epoch times corresponding to {^{t),t > 0}. Then 

I3k,n = P{T > ^k), 

and as a function of k, 13k,n is a survival function with probability mass function = {bi,n, b 2 ,n, ■■■), 
where bk,n = P{T = dk). 

Proof. We have 

n 

P{T >dk) = ^ P{T > dk\M = m)P{M = m) 

m=l 

n k 

— E ^mPC^ < fn\M = m) 

m=l 2=1 

n k 

m=l 2=1 

n 

= E - 1) = I3k,n, (7) 

m=l 

where the first equality follows from the fact that the lifetime of network is more than the arrival 
time of the kth shock if and only if the number of failed components in the time of A:th shock 
is less than m and the second equality follows because the random variable M is independent of 
Wi, W 2 , ■ ■ ■ ■ Since do = 0, and the network fails finally with one of the shocks, we have 

f3o,n = 1, lim 13k,n = 0. 

k^oo 

On the other hand, since {T > ^ {T > we get (3k+i < f3k and hence 13k,n is decreasing 

in k. Thus (3k,n, as a function of k, k = 0,1,, has properties of a discrete survival function. 
Let bn = {bi,n,b 2 ,n, ■■■) be the probability mass function corresponding to (3k,n- That is, bk,n = 
(3k-i,n — (3k,n- Then, based on ([7]), we have 

bk,n = P{T > dk-i) - P{T > dk) = P{T = dk). 

□ 
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From Proposition [H the fcth element in b^, bk,n, denotes the probability that the network fails at 
the time of occurrence of the A:th shock, i?*,. We call, throughout the paper, the vector as the 
vector of shock t-signature (ST-signature) of the network. 


In the following, we show that the reliability function of the network lifetime can be represented 
as the reliability functions of epoch times '(9*. For the counting process > 0}, it is known 

that {C{t) = k} if and only if {^k < ^ < i9fc_|_i} where 'do = 0. Using this fact, we have 

OO 

p{T>t)=Y,h,nP{m = k) 

k=0 

OO 

= X] I^KnP{^k < t < ^k+l) 

k=0 

OO 

= ^ Pk,n {P{'&k+1 >t)- P{'dk > t)) 

k=0 

OO OO 

= ^Pk-l,nP{'&k >t)- ^l3k,nP{'&k > t) 

k=l k=l 

OO 

= '^bk,nP{'&k > t). (8) 

k=l 


Remark 1. The model in ([5l), which arises in reliability theory, is known as the damage shock 
model (see [2], p. 92). Let a device be subject to shocks appearing randomly over time. Assuming 
that the device has a probability P{k) of surviving the first k shocks, k = 0,1 ,..., and N{t) denotes 
the number of shocks that the device is subject to in the interval [0,t], then the reliability of the 
device, H{t), at time t is 


"(*) = E PiN{t) = k)P{k), t > 0. 

k=0 

Various properties of this model have been explored by different authors; see, for example, m-m- 


The hazard (failure) rate of a random variable X or its distribution F with density function / is 
defined by Xf{x) = f{x)/F{x), where F = 1 — T is the survival function of X. The distribution 
function F is said to be increasing hazard rate (IHR) if F[t + x)/F{t) is decreasing in t whenever 
X > 0. From representation ([8]), the hazard rate of the network can be written as 

OO 

m = E Pk,n{t)Xk{t), 

k=l 


bk,nP{'&k > t) 

Ylj=l kj,nP{'&j > t) 
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where Afc(t) is the hazard rate ofand 

Pk,n{t) = 



It is interesting to note that Pk,n{'t) can be written as Pk(t) = P(T = 'dklT > t). This is true because 


bk^nP{^k >t) 

where the second equality follows from the fact that {'dk > t} and {T = 1 ?^} are independent. 

In the following, we make some stochastic comparisons between the performance of two networks, 
where the components of the networks are subject to failure according to different or same counting 
processes. We first use the following ordering definitions. 

Definition 2. Let X and Y be two random variables with survival functions F and G having 
density functions f and g. 


(a) X or F is said to be stochastically less than or equal to Y or G, denoted by X <st Y or 
F <st G, if F{x) < G(x) for all x. 


(b) X or F is said to be less than or equal to a random variable Y or G in hazard rate order, 

G{x) 

denoted by X <hr Y or F <hr G, if ' increases in x. 

F{x) 

(c) X or F is said to be less than or equal to a random variable Y or G in likelihood ratio order, 

q( Xi 

denoted by X <ir Y or F <i^ G, if is an increasing function of x. 

fix) 

We have now the following theorem. 


Theorem 1. Consider two networks consisting of ni and n 2 components and lifetimes Ti and T 2 , 
respectively. Suppose that the components of the ith network are subject to shocks which appear 
according to counting process > 0}, i = 1,2. Let the hn} = ^2 ■••)? * = 1,2 denote 

the ST-signature of the ith network. If CiiG ^st f, 2 it) and <st lon} then Ti <st T 2 . 


Proof Take * = 1, 2. Then, using ([5]), we have 

00 

P{T, > t) = P(?i(t) = k) 

k=0 

00 

k=0 

00 

<Y.^klPiUt) = k) = PiT2>t), 

k=0 
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where the first inequality follows from the facts that 13k ,* = 3,2 is decreasing in k and the 
assumption ^i(t) >st ^ 2 {t). The second inequality follows from the assumption that <st 


b 


( 2 ) 
n2 ■ 


□ 


Corollary 1. In Theorem^ assume that the eomponents of the two networks are subject to failure 
by shocks appear according to renewal processes {f,i(t),t > 0} and > 0}, respectively. 

Let Xij, i = 1,2, j = 1,2,..., denote the time between the (j — l)th and jth shocks in the ith 
network. Then the result of the theorem remains valid if we replace the condition f,i{t) >st 
with Xi,i <st X 2 ,i. 


Proof. Let Di^k = = 3,2,k = 1,2,.... Using Theorem 1.A.3 (b) of [T9j, the condition 

Xi,i <st X 2 ,i implies that Hi^k ^st '& 2 ,k- Now the result follows from Theorem [1] and the fact that 
for any counting process {^{t),t > 0 } with occurrence times ili,'d 2 , ■ ■ ■, we have {??„ < t} if and 
only if {f{t) > n}. □ 


Before presenting the next theorem, we give the following definition (see, m)- 

Definition 3. Let A and B be two subsets of the real line. A non-negative function K defined on 
A X B is said to be totally positive of order 2, denoted TP 2 , if for all ai < 02 , and bi < 62 , (oi S A, 
bieB,i = 1,2), 

K{a 2 , b 2 )K{ai,bi) - K{ai, 62 )^ 1 ( 02 , ^ 1 ) > 0. 


In the next theorem, we show when ST-signature vectors of two networks are hr ordered then the 
lifetimes of the networks are also ordered in hr ordering. 

Theorem 2. Assume that the assumptions of Theorem{^ are met and that the eomponents of two 
networks are subject to failure by shocks appear aecording to the same counting process {^{t),t > 0}. 
Ifhn} <hr b® and P{f,{t) = k) is TP 2 m fc € {0,1, ... } and t > 0, then Ti <hr T 2 . 

Proof Let = 3,2. The assumption b^i\^ <hr implies that 

is increasing in k. Then, according to Definition [31 it can be concluded that is TP 2 in 

k G {0, 1,...} and i € 1, 2. Thus, if P{f,{t) = k) is TP 2 in A; € {0,1, ...} and t > 0, then from basic 
decomposition formula (see, m). we get 

00 

p(T, > t) = k) 

k=l 

is TP 2 in i G {1, 2}, and t ^ 0 which in turn implies thcit T\ ^ 2 * 
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Remark 2. In Theorem [H if we assume that the components of two networks fail by shocks 
appear according to the same renewal processes based on i.i.d. r.v.s Xi,i = 1 , 2 ,..., then under 
the assumption that <hr hn^ and that Xi has increasing hazard rate, we have Ti <hr T 2 . 
This is true because when Xi has increasing hazard rate then <hr '^k+i^ k = 1,2,.. and hence, 
the required result follows from the representation ([8]) and Theorem 1.B.14 of m- Also if Xi has 
log-concave density function, then Dk <ir '&k+i,k = 1,2,.... Thus using Theorem 1.C.17 of 121, if 
bn\^ <ir and Xi has log-concave density function then Ti <iy. T 2 . 


3 A binomial based model 

In this section, we consider the shock model is presented in Section 2 and assume that the number 
of component failures at each shock follows a binomial distribution. Suppose that when a shock 
arrives each component fails with probability p. Assuming that the components fail independent 
of each other, the number of failed components in the first shock, Wi, has binomial distribution 
b[n,p), where n is the number of components in the network. Suppose that, the number of failed 
components in the ith shock, Wi, i > 2, depends only on Wi, ...,Wi-i through X]}=i 
binomial distribution b{ni,p), where n* = n — other words, assume that 


P{Wi = k)=\^ ]p^q^‘ 


k = 11,1 ,..., n 


(9) 


and for i > 2, 


2—1 


P{Wi = k\'^Wj =w) 
1=1 


^ ^ ^ k = 0, ...,n — w, w <n. 


( 10 ) 


where q = 1 — p. 

Now we can prove the following lemma. 

Lemma 2. Under the assumptions (0) and m, we have 


i=\ 


j = o,...,n. A; = 1,2,. 


Proof. We prove the lemma by induction. For k = 1 the result is true by relation ([9]). Assume that 
the result is true for k = m. That is 


P{^W,=j)=r:]{l-qn^q 

i=l 


_ rn\j rn{n—j) 
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Then, for /c = m + 1, we get 


m+1 

P{^W,=j) 


i=l 


^ P{Wm+l =j-k\Y,Wi = k)P{Y^ Wi = k) 

k=0 i=l i=l 


k=0 

3 


j — kJ \k 

k=o ~ *=0 


m—1 


„n{m+l)-j [ A / Z^i=0 g \k 

j) 

\j) 


"i+i V „(”^+i)(«-—i) 




which is the required result. 


Now, based on the model given in ([5]), the reliability of the network at time t is 

OO 

p{T>t) = Y,K,nPm = k), 

where /3 q „ = 1, and for A; = 1, 2,... 


A:=0 


n—1 


= J(i-W”-’' 

3=0 


n i—1 


E<E 

i=l j=0 


n n—1 

EEs; 

m=l j=n—m 


n\ J 


J J \n — m 


(^_iy-n+m^km 


□ 


( 11 ) 


( 12 ) 


From representation (l8|), we have 

OO 

PiT>t) = Y,bl,nP{^k>t), 
k=l 

where = (^Li,n “/^fc,n)- 

In the following, we concentrate on a special case where the shocks appear as a nonhomogeneous 
Poisson process (NHPP). Recall that a counting process > 0} is called a NHPP if the 
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survival function of arrival time 'dk of the fcth event is 

= <> 0 . *= 1 . 2 ..... 

x=0 

where A{t) = E{N(t)) = — log G{t), and G{t) is the reliability function of the time to the first event. 
The function A(t) is called the mean value function (m.v.f.). For more details on the properties of 
NHPP and related processes, one can see, for example, m- 

Let us look at the following example. 


Example 2. Consider a series network consisting of n components. Suppose that the network is 
subject to shocks which appear according to a NHPP with m.v.f. A(t) = — logG(t). Then under 
model (fTTT) and noting that the t-signature of a series network is = (1,0,0,... , 0), we can easily 
see that 

= A: = 0,1,2,.... 

Hence, the reliability of series network is given by 


PiT>t)='£Q""pm = k) 

k=0 


k=0 


(Mt))” 

k\ 


=e ^ 
k=0 


(g"A(t))'= 

kl 


=(G(t)y-p 


Note that if n, the number of components of the network, gets large then the reliability of the 
network tends to G(t). 


In the sequel, we explore some aging properties of the network lifetime. First, recall that a 
distribution F is said to be increasing hazard rate average (IHRA) if is decreasing in 

t > 0. It is well known that the IHR property implies the IHRA (see [2]). 

We have the following lemma. 

Lemma 3. ^ is IHRA. 

Proof. In order to prove the result, we must show (/3^ * is decreasing in k for A: = 1, 2,... . Note 
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that /3^ ^ can be rewritten as 


n « 

j-i 


1 - uY-^ 

k B{j,n- j + l) 


du, 


(13) 


which is clearly an increasing function of q^. It is clear from (|1‘2I) that ^ is a static reliability 
function of a network. If we write /3^ „ = h{q^), where h is the reliability function of the network, 
then by choosing a = in Theorem 2.5 of Section 4 of [2], we conclude that 




which is equivalent to say that 


In)-^ > 


t+l,n. 


1 

I k + l 


This completes the proof of the lemma. 


□ 


The following example shows that, although /3^ ^ is always IHRA, but it is not necessarily IHR. 

Example 3. Consider a bridge network pictured in Figure [5J It can be seen that the t-signature 
of this network is as = (0, 0). In order to show that ^ is IHR we have to show, 

based on the definition of IHR distributions, that is decreasing in k. 


®\ ^ )® 


Figure 2: The bridge network. 

Figure [3] shows the plot of for this network where q = 0.5. As the plot shows, this ratio is 

k'k 

not decreasing for all values of k, hence is not IHR. 

Theorem 4.1 of m implies that if P{^{t) = k) is TP 2 in t € (0,oo), and k G {0,1,...} and 
(E(a«W))i is decreasing in t for a G (0,1), then based on the fact that ^ is IHRA we get that T 
is also IHRA. The following theorem shows that under the condition that ^{t) is NHPP, the network 
lifetime T is IHRA if the distribution function of the arrival time of the first shock is IHRA. 

Theorem 3. Consider a network consisting of n eomponents with lifetime T. Suppose that the 
components of the network is subjeet to failure by shocks that appear aceording to a NHPP with 
m.v.f. A{t) = — logG'(t). If G is IHRA, then T is IHRA. 
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Figure 3: The plot of for the bridge network. 


Proof. For ki < k 2 , 


= k2) _ ^ 

pm = ki) k2r 


is increasing in t and hence = k) is TP 2 in k and t. On the other hand, for a € (0,1) 

Z—/ ri' 


n=0 


ni 


=e-Mt) ^ ^ g-A(t)(l-a) 


n=0 
, 1—a 


ni 


= {G{t)y 

If G{t) is IHRA, then (G(t)) * is decreasing in t and hence 

is decreasing in t. Hence the result follows from Theorem 4.1 of Gottlieb |10j . 


□ 


In the next theorem the stochastic relationships between t-signature vectors and the lifetimes of 
two networks are investigated. 

Theorem 4. Consider two networks with lifetimes Ti and T 2 and t-signature vectors s[ = (sj ..., sj „) 
and S 2 = respectively. Suppose that the components of the ith network is subject to 

failure by shocks appear according to NHPP with m.v.f. Aj(t) = — logGi(t), i = 1,2. Assume that, 
upon arriving the shocks, the components of the ith network fail with probability pi, i = 1,2. 
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(a) If Pi > p2, Gi <st G2 and s[ <«* then Ti <st T2. 

(b) If Pi = p2, Gi =st G2 and s[ <hr then Ti <hr T2. 

Proof. Let = Yl^k+i where is the jth element of ST-signature associated to the ith 
network and qi = 1 — pi, i = 1, 2. 


(a) Let {f,iit),t > 0} be the NHPP with m.v.f. Aj(t) = — logGj(t), i = 1,2. Supose that 
9j,niQi) = iLgk du, i = l,2. Using ([I3D, it can be seen that ), 

where gj^nili) is an increasing function of gf. Hence 


i=i 

n 

n 


^ '^^lj9j,n{92) = Pi, 
i=i 


*( 2 ) 

n 


in which the first inequality follows from the fact that sj <st S 2 and gj,n{Q^) is increasing in 
j and second equality follows from the assumption pi > p 2 which implies gj^niQi) < 9 j,n{g 2 )- 
Also, Gi <st G 2 implies ^i(t) >st ^ 2 (i)- Then the result follows from Theorem [TJ 

(b) It is easy to see that (j)(l — is TP 2 in k and j. Also, s[ <hr sj implies that Sfj 

is TP 2 in i and j. Therefore from basic decomposition formula (see m), 

n-l , ^ 

j=o ^ 


is TP 2 in A: € {0,1,...} and i G {1,2} which implies <hr The proof is complete 

based on Theorem [2l 


□ 

Example 4. Consider again Example [3j Let the network be subject to shocks that appear accord¬ 
ing to a NHPP with m.v.f. Aft) = —logGft) and in each shock, each link fails with probability 
0.1. We are interested in assessing the reliability of the network in the cases where the time to the 
first shock has either an exponential distribution with a constant hazard rate of 1 (Exp(l)), or a 
Weibull distribution with shape parameter 2 and scale parameter 1 (W(2,l)), or a linear hazard 
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distribution (L(l,l/2)). The survival functions of these distributions, respectively, are given as 

G'i(t) =exp(—t), t > 0, 

G 2 (t) =exp(-t^), t > 0, 

G 3 {t) =ex.p{—t — t^), t > 0. 

It can be easily shown that L(l,l/2) is stochastically less than both Exp(l) and W(2,l). Hence, 



Figure 4: The plot of network reliability in Example IH 

as Figure 0] reveals, based on Theorem 01 the reliability of the network for the L(l,l/2) case is 
less than that of the cases of Exp(l) or W(2,l). It can be easily seen that Exp(l) and W(2,l) are 
not stochastically ordered. Also, the plot shows that the network lifetimes are not stochastically 
ordered. 

4 Network reliability under fatal shocks 

In this section, we assume that each shock is fatal for the network. That is, when a shock arrives it 
leads to failure of at least one component. Let fatal shocks occur according to a counting process, 
> 0}, at random time instants Qi, g 2 , ■ ■ ■ ■ It is clear that the network finally fails by one of 
the fatal shocks. In order to obtain the reliability function of the network, in such a situation, first 
we obtain P{T = Qi), i = 1,... n. Consider a network consists of n components. It can be shown 
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that the number of ways showing the order of component failures is n* given in Lemma [H Then, 
under the assumption that all ways of the order of component failures are equally likely, we have 

s* = P{T = Qi) = ^, i = 
n* 

where is the number of ways of the order of component failures in which ith fatal shock causes the 
network fails. It is obvious that s* just depends on the structure of the network. In the following 
example, we compute s* = (sj^,..., s*). 

Example 5. Consider again Example [T] Let vr denote the order of link failures in the network and 
r(7r) the shock number that caused the failure of the network. All possible vr and corresponding 
r(7r) have been presented in Table [2j It is clear that 

st = PiT = Qi) = ^ s*^ = P{T = Q2) = ^, S3 = P{T = Q^) = 0. 

That is, s* = (^, ^, 0). 


Table 2: All possible vr and corresponding r(7r) 


TT 

r(7r) 

TT 

r(7r) 

TT 

r(7r) 

(1,2,3) 

1 

({1,3},2) 

1 

({1,2,3}) 

1 

(1,3,2) 

1 

({2,3},1) 

1 



(2,1,3) 

2 

({1,2},3) 

1 



(2,3,1) 

2 

(3,{1,2}) 

2 



(3,1,2) 

2 

(2,{1,3}) 

2 



(3,2,1) 

2 

(1,{2,3}) 

1 




From the fact that s* does not depend on the random mechanism of the component failures, we 
obtain the reliability function of the network as 

n 

PiT > t) = ^ P(r > t\T = Qi)PiT = g,) 

i=l 

n 

= '^s*P{gi > t\T = gi) 

i=l 

n 

= ^s*Pigi>t). (14) 

i=l 

From the fact that C(i) = k if and only if t it can be seen that 

n—1 

PiT>t)=^s:pm = i) ( 15 ) 

1=0 
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where S* = E]=i+i s*j- 

Remark 3. It is noted that the representations (|14p and ()15p are similar to representations ([SP 
and m, respectively. Hence, the results obtained based on ([8|) and (l5|) in Section [2] are valid for 
the fatal shock model. 
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